AN INVERSE PROBLEM IN 
QUANTUM STATISTICAL PHYSICS 



FLORIAN MEHATS AND OLIVIER PINAUD 

Abstract. We address the following inverse problem in quantum statistical 
physics: does the quantum free energy (von Neumann entropy + kinetic energy) 
admit a unique minimizer among the density operators having a given local 
density n{x)? We give a positive answer to that question, in dimension one. 
This enables to define rigourously the notion of local quantum equilibrium, or 
quantum Maxwellian, which is at the basis of recently derived quantum hydro- 
dynamic models and quantum drift-diffusion models. We also characterize this 
unique minimizer, which takes the form of a global thermodynamic equilibrium 
(canonical ensemble) with a quantum chemical potential. 



1. Introduction 

We deal with a question which is at the core of recently derived quantum hy- 
drodynamic models based on an entropy minimization principle |13| Let a 

given density of particles n{x) > 0, can we find a minimizer of the quantum free 
energy among the density operators g having n{x) as local density, i.e. satisfying 
the constraint p{x,x) = n{x), where p{x,y) denotes the integral kernel of g? 

This question arises in the moment closure strategy initially introduced by De- 
gond and Ringhofer in [13j in order to derive quantum hydrodynamic models from 
first principles. Let us briefiy review this theory (for more details, one can refer to 
the reviews |12l llOj). The quest of macroscopic quantum models is motivated by 
applications such as nanoelectronics, where affordable numerical simulations of the 
electronic transport are necessary while the miniaturization of devices now imposes 
to take into account quantum mechanical effects in the models, resulting in a higher 
simulation cost. At the microscopic level of description, the Schrodinger equation 
and the quantum Liouville equation are numerically too expensive, which motivates 
the derivation of models at a more macroscopic level. In the classical setting, the 
relationships between microscopic (kinetic) and macroscopic (fluid) levels of descrip- 
tion are fairly well understood by means of asymptoti c analysis, see for instance 
[18^ [171 ES]- In particular, it is known that the understanding of the structure of 
the fluid model relies on the properties of the collision operator at the underlying 
kinetic level. Indeed, collisions are the source of entropy dissipation, which induces 
the relaxation of the system towards local thermodynamical equilibria. The free 
parameters of these local equilibria are the moments of the system (e.g. local den- 
sity, momentum and energy) and are driven by the fluid equations. Arguing that 
the derivation of precise quantum collision operators is a very difficult task, while 
only the macroscopic properties of such operators is needed in our context, Degond 
and Ringhofer have grounded their theory on a notion of quantum local equilibria. 
To do so, they have generalized Levermore's moment approach |27| to the quantum 
setting. The idea consists in closing the system of moment equations by defining a 
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local equilibrium as th e minimizer of an entropy functional (say, the von Neumann 
entropy) under moment constraints. 

In this approach was adapted so as to describe systems in strong inter- 

action with their surrounding media and obtain quantum macroscopic models by 
applying a diffusive asymptotics. The most simple of these models, the quantum 
drift-diffusion model, was studied numerically in [Ml [7] and the simulation results for 
one-dimensional devices such as resonant tunneling diodes were encouraging. This 
model is based on the most elementary constrained entropy minimization problem. 
Indeed, in this case, the local quantum equilibrium at a given temperature, also 
called quantum Maxwellian, is defined as the minimizer of the quantum free energy 
subject to a local constraint of prescribed density. Note that not only the total 
number of particles is fixed, as in the usual quantum statistics theory (for the so- 
called canonical ensemble), but also the local density n(x) is imposed at any point x 
of the physical space. This problem has been studied formally in citeQET and the 
Lagrange multipliers theory lead to the existence of a quantum chemical potential 
A{x) such that the solution of the minimization problem is a density operator of 
the form 

Q = exp[ . (1.1) 



Remark that the difficulty in this problem lies in the fact that its solution will 
depend on its data in a global way. The similar problem in classical physics, i.e. 

I 1 2 

reconstructing f{x,v) = exp( — (^^^ — \- A{x))) from its density n{x) = f f{x,v)dv, 
is very simple and the chemical potential, given by A{x) = — logn(x) + | log(27r), 
depends on n{x) in a local way. Here, due to the operator formalism of quan- 
tum mechanics, which is not commutative, the density and the associated chemical 
potential are linked together by a non-explicit formula, and in a global manner. 

To end this short presentation, let us also recall that this quantum drift-diffusion 
model displays formally several interesting properties: it dissipates a quantum fluid 
entropy, which indicates that it should be mathematically well-posed, and it can be 
related to other known models after some approximations (for instance, semiclassi- 
cal expansions on the quantum drift-diffusion system enable to derive the density- 
gradient model). Besides, a whole family of quantum fluid models were derived 
by several authors, based on the same entropy minimization principle: quantum 
Spherical Harmonic Expansion (QSHE) models [1], quantum isothermal Euler sys- 
tems [22] (8j, quantum hydrodynamics \23\ [9], models with viscosity [Gj \2T\ 120] . 
quantum models for systems such as subbands [32] or spins [3]. Nevertheless, one 
has to put the emphasis on the fact that all these studies rema in yet at a formal 
level. Even the notion of local quantum equilibrium has only been deflned formally 
and this problem of entropy minimization under local constraints is widely open. 

The aim of this paper is to make a flrst step towards the rigorous justification of 
these models, by studying the quantum entropy minimization principle in the most 
simple situation, in the case of a density constraint. We work in dimension one, in 
a finite box with periodic boundary conditions. Our main result, Theorem 12. H is 
presented after a few notations in the next section. We show that, in an appropriate 
functional framework, the quantum Maxwellian is properly defined, i.e. that to any 
density n(x) > corresponds a unique density matrix g minimizing the free energy. 
Moreover, we prove that g actually takes the form (jl.ip . where A(x) is the quantum 
chemical potential (in the sequel of the paper, the temperature T will be set to 1). 
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Let us now make a remark. One can see on the formula (jl.ip that the quantum 
Maxwehian reads as the global equilibrium canonical ensemble associated to the 
Hamiltonian —A + A(x), where the chemical potential A(x) is seen as an applied 
potential. Hence, our problem can be reformulated as the following inverse problem 
in quantum statistical mechanics. Let a system at thermal equilibrium with a sur- 
rounding media at a given temperature, in a certain potential. Can we reconstruct 
the potential from the measurement of the density at any point? This problem has 
been much less studied than more standard inverse problems such as the inverse 
scattering theory (reconstructing the potential from its scattering effects) or the 
inverse spectral problem (reconstructing the potential from the spectrum of the as- 
sociated Hamiltonian). Nevertheless, one can quote at least two references where 
similar inverse problems have been investigated. In [26j (see also the series of ref. 5 
therein), a practical method for reconstructing potentials from measurements was 
settled using Feynman path integrals and, in [19j, a close problem for quantum spin 
systems was studied. 

The outline of this paper is as follows. In Section [21 we define the functional 
framework of the paper and state our main theorem. In Section [3l we study the 
entropy and the free energy and give some useful results for the sequel. In Section 
m we prove the existence and uniqueness of the minimizer Q[n] associated to a 
density n. SectionOis devoted to the characterization of Q[n] via the Euler-Lagrange 
equation for the minimization problem. To deal with the constraint, we introduce 
a penalized problem. 

Future developments of this work will involve several directions. A first exten- 
sion will concern the investigation of other spatial configurations: other boundary 
conditions, whole-space case, or space dimension greater than one. We will also 
investigate the entropy minimization problem with constraints of higher order mo- 
ments. As in the case of classical physics, it might lead to ill-posed problems and to 
delicate problems of moment realizability |24| . Another interesting question which 
remains to be solved concerns the quantum evolution: can we define an evolution 
for a quantum Liouville equation with a BGK-relaxation operator based on the 
local equilibria defined in this paper, as for instance in [Ij for other relaxation op- 
erators? This issue is linked to the possibility of rigourously deriving the quantum 
drift-diffusion model. 



2. Notations and main result 

Let us describe the functional framework of this paper. The physical space that 
we consider is monodimensional and bounded. The particles are supposed to be 
confined in the torus [0, 1], i.e. with periodic boundary conditions. We consider the 
Hamiltonian 

H = -^ 

on the space L^(0, 1) of complex- valued functions, equipped with the domain 

D{H) = 1^ G R\0, 1) : ^(0) = n(l), ^(0) = ^(1)| . 

The domain of the associated quadratic form is 

H^,, = {nGHi(0,l): u{0) = u{l)} . 
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Its dual space will be denoted Hpg^,. Remark that one has the following identification: 



Vn, u G Hp 



per) 



du 



dx 



(2.1) 



L2 



We shall denote by Ji the space of trace class operators on L^(0, 1) |30| [33] and by 
J2 the space of Hilbert-Schmidt operators on L^(0, 1), which are both ideals of the 
space >C(L^(0, 1)) of bounded operators on L^(0, 1). We denote by /C the space of 
compact operators on L^(0, 1). 

A density operator is defined as a nonnegative trace class self- adjoint operator 
on L^(0, 1). Let us define the following space: 

£ = [q(^ Ji, Q= Q* and VH\Q\^fH G Ji} . 
This is a Banach space endowed with the norm 

ll^ll^- = Tr + Ty{Vh\q\Vh). 
For any Q ^ £, the associated density n[^] is formally defined by 

n[Q]{x) = p{x,x), 
where p is the integral kernel of g satisfying 



£L\0,1), q{<P){x) 



pix,y)4>iy)dy. 



The density n[g] can be in fact identified by the following weak formulation: 

V$GL~(0,1), Tr($^)=/" ^{x)n[0]{x)dx, (2.2) 

Jo 

where, in the left-hand side, $ denotes the multiplication operator by ^, which 
belongs to £(L^(0, 1)). If the spectral decomposition of g is written 

00 



k=l 



then we have 



Tv ^\g\^ = \\y/Hy^\\\% =^ 



Pk\ 



k=l 



dx 



2 

L2 



n 



\n 



MIIli <5^|/0fc| =Tr|i?|. 



(2.3) 



(2.4) 



k=l 



k=l 



Moreover, by the Cauchy-Schwarz inequality, n[g] belongs to VF"'^'^(0, 1) with peri- 
odic boundary conditions, and we have 



dx 



1/2 



< 2\\gfl^ (TTy/H\g\y/HY'^ < C\\g\\£. 



The energy space will be the following closed convex subspace of £: 

£+ = {ge£ : g>0}. 
On £^ we define the following free energy: 

F{g) = Tr(^log(^) - g) + Tr{y/H g^/H). (2.5) 
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We will see in Section [3] that F is well-defined and continuous on £j^. If G 
then the Cauchy-Schwarz inequality applied to (|2.4p gives 



d_ 

dx 



< 



\k=l 



dx 



1/2 



Hence we have n[Q] G H^g^ and, using (|2.3p . we get 



L2 



1/2 



(2.6) 



Recall also the following logarithmic Sobolev inequality for systems, proved in |14) 
and adapted to bounded domains in [15]: for all g £ £^ we have 



Tr Q log + Tt{Vh gy/H) > f n[g] log n[g]dx -\ ^ — ^ Tr g. 

Jo 2 



(2.7) 



This inequality, coupled to (|2.6p which gives logn[£»] G L^{0,1), implies that 
Tt glog g is bounded for all G 

Our main result is stated in the following theorem. 

Theorem 2.1. Consider a density n G H^g^, such that n > on [0,1]. Then the 
following minimization problem with constraint: 



mini<'(^) for g £ such that n[g\ 



(2.8) 



where F is defined by (|2.5p . is attained for a unique density operator g[n], which has 
the following characterization. We have 

g[n]=exp{-{H + A)), (2.9) 

where A belongs to the dual space Hp^,. of H^g^, and the operator H + A is taken in 
the sense of the associated quadratic form 

2 



dip 



dx 



L2 



+ {A,\ip\%-i 

/ riper .ripg 



(2.10) 



From (j2.9p . it is possible to obtain a formula for A. Such formula is given in 
(j5.50p . The following remark shows that the functional space H~g^, for the quantum 
chemical potential A{x) is optimal. 

Remark 2.2. For a given A G H~g^, we prove further -see subsection 15.21 Step 4 of 
the proof- that the operator H+A{x) (in the sense of quadratic forms) is self-adjoint 
and has a compact resolvent. Moreover, the associated quadratic form is a form- 
bounded perturbation of u i— )• ||u'||^2) so that H + A{x) can be diagonalized on 
and its A;-th eigenvalue has an asymptotic behaviour of the form C/c^. Therefore, 
the norm of the associated eigenvector (pk is bounded hy Ck. Consider now the 
operator g = exp{—{H + A)) and the associated density n{x). By using the decay 
of the exponential, one can see that the series in (|2.4p is converging in H^. If we 
assume that A belongs to the Sobolev space H'^ , where — 1 < s < 0, then by elliptic 
regularity one has 4>k G H^"'"^, and the series (|2.4p will converge in this Sobolev 
space, so we deduce that n G H*"'"^. Hence, if n belongs to H^ but does no t belong 
to any H*, s > 1, then we have A G H~g^ and A cannot be more regular. 
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Remark 2.3. The two main limitations of this theorem, the strict positivity of n 
and the one-dimensional setting, are not essential for the first part of the theorem, 
the existence and uniqueness of the minimizer. This first result will be extended 
in a forthcoming work. However, these assumptions are essential in our proof of 
the second part of the theorem, the characterization of the minimizer. Indeed, the 
strict positivity of n is crucial in Subsection 15.21 Step 3, see e.g. Eq. (|2.9p and 
the argument after ()5.5ip . Moreover, the one-dimensional framework implies by 
Sobolev embeddings that Hpg,. is a Banach algebra, which enables to define the 
above quadratic form Qa in (|2.10p . 

Outline of the proof of Theorem 12.11 The existence of the minimizer Q[n] of the 
constrained problem (12. 8p is obtained by proving that minimizing sequences are 
compact and that the functional is lower semicontinuous. Compactness stems from 
uniform estimates that enable to apply Lemma 13. H whereas the lower semiconti- 
nuity comes from (|3.2p in Lemma 13.11 and from Lemma 13.31 The uniqueness of the 
minimizer is a consequence of the strict convexity of the entropy (proved in Lemma 

mi- 

In order to characterize the minimizer g[n] of (|2.8p . we need to write the Euler- 
Lagrange equation for this minimization problem. This task is difficult because the 
constraint n[Q] = n is not easy to handle when perturbing a density operator. We 
circumvent this difficulty by defining a new minimization problem with penalization, 
whose minimizer will converge to Q[n]. Next, the Euler-Lagrange equation for 
the penalized problem reads 

V& (log(&) + H + A,) 

From this equation, in order to prove that 

& =ex.p{-{H + Ae)) , 

we show two important intermediate results, relying on the fact that is a min- 
imizer and on properties of the von Neumann entropy Tr(£)log£) — g). First, we 
prove that the kernel of g^ is reduced to {0}. Second, we prove that the family {(pp) 
of eigenfunctions of g^, which is a Hilbert basis of (0,1), is in fact dense in Hi„^. 
This enables to prove that {(pp) is the complete family of eigenfunctions oi H + A^, 
and to identify the associated eigenvalues. Finally, using the two assumptions dis- 
cussed in Remark 12. 3| we are able to prove that A^ converges in the Hpg^ strong 
topology, which is sufficient to pass to the limit as the penalization parameter e 
goes to zero, and conclude the proof. 

3. Basic properties of the energy space and the entropy 

In this section, we prove a few basic results on the energy space £"+ that will be 
used in the paper. 

Lemma 3.1. Let {gk)k£N be a bounded sequence of£^. Then, up to an extraction 
of a subsequence, there exists g £ £^ such that 

Qk ^ Q in J\ and y/g^ —5- in J2 as k ^ +00 (3.1) 

and 

Tr{VHgVH) < liminf Tr(\/FaV^). (3.2) 
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Furthermore, if one has 

TT{^/Hg\^) = lim Tv{^/Hgk\^) 

A,'— >+oo 

then one can conclude in addition that 

V^-y/^ —5- ^J~H.^fg}^ in Ji as k ^ +00. (3.3) 

Proof. Step 1: weak-* convergence in J\. First notice that the boundedness of g^^ in 
£,j^ imphes by (|2.3p that the operator ^/H^fgk is bounded in the Hilbert space J^2- 
Moreover, since g^ and \fHg^\fH are positive and bounded in we can extract 
subsequences such that and ^/Hg^-jH converge in the J\ weak-* topology, that 
is, there exists two positive trace class operators g, A such that, for all compact 
operator K £ IC, 

Tr(Ka) ^ TT{Kg) ; TrixVHgkVH) Ti{KA). 
By application of Proposition 3.12 of [5], we have 

\\A\\j^ < liminf \\^/Hgk\^\\j^ = liminf Tr(\/FaV^)- (3-4) 

Step 2: Identification of A. We show that A = \/Hg\fH. Indeed, let K = {^/H + 
K' {^/H + with K' compact. Using the cyclicity of trace with the bounded 
operators yJ^\fH and \fH^fg}^^ we get 

TriKVHgk^) = T,{^kVH K Vh ^t). 

= Tt{^K'^) - Tt{^K'{^ + I)-^^) 

+ Tii^iVW + I)-^K'{^/H + 1)-^^), 

= Ti{K'gk) - Tt{K'{^/H + I)-^gk) - Tt{{Vh + ly^K'gk) 

+ + I)-^K'{Vh + ly^gk)- 

Since gk ^ Q the J^i weak-* topology, and since (\/^ + K' {^/H + , 
K'{^/H + I)^^ , {\fH + I)^^K' are compact operators, we have 

TI{K^/HgkVH) Ti{K' g) - Tt{K'{Vh + I)-^ g) - Tr{{VH + ly^K'g) 

+ Tii{^/H + I)-'K'{Vh + I)-^g) 
= TTiK'iVH + I)-^A{VH + I)-^). 

We thus obtain 

+ = g-{^/H + I)-^g-g{VH + I)-^ 

+{^fH + I)-^ g{VH + 
and it follows that A = \fHp\fH. In particular, (|3.4p yields (|3.2p . 



Step 3: weak convergence in J\. Let us prove now that g^^ converges weakly in 
that is, for all bounded operator a G £(L^(0, 1)), 

Tr(o-£ifc) Tr(cj^). 
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We have 

Tr((T£»fc) = Tr(^/^cj^/^) 

= Tr(^(^/F + + I)-^<j{Vh + I)-\Vh + /)^) 

= Tt{{Vh + iy^a{VH + I)-\^/H + I)Qk{VH + I)) 

= Tr((\/1^ + I)-^a{^/H + ly^VHgkVH) + Tr((\/¥ + 1)-^^^) 

+ Tr(a(^/:H^ + /)-i£>fc) - Tr((^/F + I)-^a{VH + 

Since {^/H + 1)^^ a{-\fH + 1)^^ , {^/H + I)^^a and cr(-v/^+-^)^^ are compact, since 
Qk ^ Q and y/HQk\/H — s- \/Hq^/H in ^Ti weakly-*, we can pass to the hmit in the 
latter expression and obtain the weak convergence of Qk- 

Step 4: strong convergence in J\. To obtain the strong convergence in it suffices 
now to apply a result of |33) that we recall here (specified to the case of ) . 

Theorem 3.2 (Theorem 2.21 and addendum H of [33])- Suppose that A^. ^ A 
weakly in the sense of operators and that \\Ak\\j-^ — )• HAHj-^. Then \\Ak — A\\j-^ — )• 0. 

One can indeed apply this result since converges weakly in to g (which 
implies the weak operator convergence), with convergence of the respective norms: 

WekWji = TTgk^TTg= \\g\\j^. 

This implies that the convergence of gk in J7i is strong and the first part of the 
lemma is proved. 

Step 5: strong convergence of in J^- We have 

\Qk - q\c{J?) < \\Qk - qWji- 

Moreover, it is known that the norm convergence oi g^ > to g > implies the 
norm convergence of to (see e.g. [30j). We claim that in fact we have 

\/^ in J2- (3.5) 

To prove this fact, since J'2 is a Hilbert space and since 

WVSkWji = Tra Tr£) = \\g\\%, 

it suffices to prove that ^ in weak. Let cr £ ^2- One can choose a 
regularizing sequence with finite rank such that 

ak ^ cr in ^2 as n — )• +00. (3.6) 

For all n,m G N, we have 

I Tr(,/^- ^)(t| < I Tr(^/^- ^)fTr„| + I Tr(^/^ - ^)(fT - 

< WVOk - ^/Q\\c{L2)\\crm\\J^ + {\\y/gk\\j2 + II J'2 ) Ikm - (^Wji 

< \/^ll£(L2)||o-m||j-i + C\\am - cr\\j^, 

which implies by ()3.6p that Tr(y^ — y/g)cr — as n — )■ +00. This means that 
y^ ^ y^ in J'2 weak, which ffiially implies (|3.5p . 

Step 6: strong convergence of y/HyJ^ in J'2- From now on, we assume that one 
has in addition the following convergence: 

W^H^fj^ = Tri^gkVH) ^ Tv{VHg^) = \\^^\\%. 
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Consequently, if we prove that 

^/Hy/Ok ^J~H^\fQ in J2 weak, 

this weak convergence in the Hilbert space J2 will be in fact a strong convergence. 
To this aim, we consider a ^ J2 and, for all e S (0, 1), we decompose 

Tiia^^k) = Tr(a(l + e^y^VW^) + Ti{a{l - (1 + eVH)~^)^^) 

= TT{aVH^) + Tr{aeVH{^ - ^)) + Tr((a, - a)VH^) + Tr((a - a,)^/Hy^) 

with (Tg = a{l + £\fH)^^ . For all e > 0, the operator (1 + £y/H)^^^/H is bounded 
on L^(0, 1), so as^/H belongs to ^2 and the previous step implies that 



lim Ti{ae^H{^- ^)) = Q. 

k—^+oo 

Now we write 

Tiiia,-a)VH^) + Tri{a-as)VH^)\ < \\ct - a^Wj, + \\Vh^\\j,^ 

< C\\a-ae\\j2. 

Therefore, if we prove that ag converges to a in J'2 as e ^ 0, we will have 

TriaVHy/gj^) ^Tr{aVH^) as A; ^ +00 

and the proof of the lemma will be complete. 

Introduce the eigenfunctions and eigenvalues (ep,/ip)pgN* of the operator H, 
which has a compact resolvent. We have 



p=i 



which converges to ||c||^2 ^ ^ ~^ ^ comparison theorem. Similarly, for all 
If e 1/^(0, 1), we deduce from the convergence in of the series ^ fpCp, where 
ipp = Jq ip{x)ep{x)dx, that 

00 ^ 

(1 + eVHY^ip = ip„ep f m. L^(0, 1) as e 0. 

^ 1 + e^/lu, 



p=i - ' - Vi^p 

Therefore, cjg and a* = {1 + e^/Wy^a converge strongly to u as e — )• and one 
can apply Griimm's convergence theorem (see |33J , Theorem 2.19), which proves 
the convergence of cjg to a in J'2. The proof of Lemma [3.11 is complete. □ 

Lemma 3.3. The application g 1— )• Tr(£)log g— g) possesses the following properties, 
(i) There exists a constant C > such that, for all g £ we have 

TT{g\ogg- g)>-C(TrVHg^y^'^ . (3.7) 



(a) Let gk he a hounded sequence of £^ such that g^ converges to g in Ji, then 

gk log gk — gk converges to g log g — g in Ji. 

(Hi) The application g 1— t- Tr(£)log g — g) is strictly convex on 
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Proof. Step 1: proof of the inequality ()3.7p . We shall use the following inequality, 
deduced from Lemma [A . 1 1 which is proven in the Appendix: there exists C > such 
that, for all q ^ £j^, 

^p^XplQ] < CTtVHqVh, (3.8) 
p>i 

where we have denoted by (Ap[£)])p>i the nonincreasing sequence of nonzero eigen- 
values of Q (this sequence is finite or infinite). The function s i— )• /3(s) = slogs — s 
is negative on [0, e] and positive increasing on [e, +oo). Let 

[slogs - s\ 
Ci = sup < +00. 

se[0,e] 

Let ^ G <?+ and denote by (Ap[f)])p>pp the eigenvalues of g that belong to the interval 
(0, e]. We have 



-Tr/3(^>) < ^ < E V^pis] 

p>po P>PO 

1/2 / \ 1/2 



< cJ^p^xm) (Ei 

\p>Po / \p>Po , 



C / 



1/2 



which proves ()3.7p . 

Step 2: proof of (ii). Consider a sequence qi^ bounded in such that g in 

Ji. Let M = sup^ l|f?fc||£{L2) < +00. There exists a constant Cm > such that 

Vs G [0,M], Islogs - s| < Ca/s^/^ 
Thus, for all e > 0, denoting again /3(s) = slogs — s, we get 

E \/^i^piek])\ < Cm Y1 iMQk]f^<CMe'/^ E i^piSk])'^' 

Ap[£>fe]<£ Ap[efc]<£ Ap[efc]<e 

1/2 

< cm.^/MEp%[^^]) (Ei 

Vp>i / Vp>i ^ 

< Ce'/^ ( Tr y^Qk^) < Ce'/\ 



where C is independent of k and where we used ()3.8p . The same inequality holds 
for the limit g. Let e > and let us decompose 

/3(s) = /3i(s) +/32(s) = (/31.<,)(s) + (/31s>.)(s). 

For all e, one has 

TT\p{gk)-m\ < Ti\p,{gk)\+TT\^i{g)\+Ti\p2{gk)-He)\ 

< Ce^/^ + Tv\(32{gk)-P2{Q)\, 

so that the result is proved if we show that fi2{Qk) converges to (32{q) strongly in 
Ji. According to Theorem 13.21 it is enough to prove that (32{Qk) converges weakly 
to I32{q) in -Ji and that ||/32(^>fc)|| ||/?2(^?)|| J^i to obtain the strong convergence 



11 



in J\. We prove first the weak convergence. To this aim, we choose e such that 
7^ £ for ah p G N* and denote 

= max {p : Xplg] > e} . 

According to Lemma [A. 21 we have 

Xp[Qk] ^ Xpie], Vp>l, (3.9) 

and we can choose k large enough so that we have 

-^pi^fe] > ^ foi' aU p < A^ and Xp[Qk] < £ for aU p > N. 

Besides, foUowing again Lemma [A. 21 we choose some eigenbasis {4'p)peN* and (</>p)pgN* 
of and Q, respectively, such that 

VpGN*, lim 11,/.^ - </.p||i2 = 0. (3.10) 

Then, the actions of f32{Qk) and l32{g) on any tp £ L^(0, 1) read 

N N 

/32igk)v = Y,(^(^pi0k]){(l)p,^)(l)p ; (32ig)ip = ^/3{\p[Q]){(j)p,ip)(pp, 
p=i p=i 

where (•,•) denotes the L^(0, 1) scalar product (taken linear with respect to its 
second variable and anti- linear with respect to its first variable). Therefore, for any 
bounded operator B, 

N N 

TT{(32{0k)B) = Y,(^iK[ek]){4>p,B4>'p) ^ ^/3(Ap[^])(0p,i?,^p) = Tr{(32{Q)B) , 
p=i p=i 

thanks to (|3.9p . (|3.10p and the continuity of the function This proves the weak 
convergence of /32(i?fc) in J^i- Regarding the convergence of the norm, we have 
directly 

N N 

mgk)U=Y.\f^{Xp[gk])\^Y.\^{Xp[Q])\ = mg)\\j,, 

p=l p=l 

and item (ii) is proved. 

Step 3: proof of the strict convexity (Hi). We recall first the Peierls inequality [33) : 
let {ui)i>i be an orthonormal basis of 1/^(0, 1), whose scalar product is denoted by 
(•, •); then, setting /3(s) = slogs — s, we have 



j;/3((n,,^^x,)) <Tr(/3(£>)). (3.11) 



i>l 



Indeed, denoting by {\i,4)i)i>i the spectral elements oi g £ we have 

Since X]j>i |(0ji^^j)P = Ij it follows from the Jensen inequality that 

/3((tx„^^/,)) = /3 K;A,-|(0„ni)n <^/3(A,-)|(0i,ni)|2. 
\j>i / i>i 

Summing up the latter relation with respect to i and using the relation X^j>;^ | {(j)j,Ui 
1, the Peierls inequality (|3.1ip follows. 
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Consider now qi, Q2 in £"+ such that qi ^ Q2. Let t S (0,1) and denote by 
the spectral elements of the operator tgi + (1 — t)Q2- Then 

Tr (/3(t^i + (1 - t)e2)) = = E + (1 - t)Q2m). 

i>l i>l 

There exists at least one index such that {i^io , giipig) 7^ ii^io ^ Q2ipio) ■ Indeed, if 
not, we would have gi = g2 since {■ipi)i(zfq* is an orthonormal basis of L^{0, 1). Since 
/3 is strictly convex, it thus comes, 

Sl A) + (1 - t)(,^^, g2 A)) < Yl ^i)) + (1 - Q2 ^i))] . 

Using the Peierls inequality ()3.1ip to control the right hand side, it comes finally 

TV {P{tgi + (1 - t)g2)) < t Tr {(3{gi)) + (1 - t) Tr {/3{g2)) , 
which yields the strict convexity of the functional. □ 



4. Existence and uniqueness of the minimizer 

In this section, we prove the first part of our main Theorem 12.11 More precisely, 
we prove the following proposition. 

Proposition 4.1. Consider a density n{x) such that n > on [0, 1] and n G H^g,,. 

Then the minimization problem with constraint 

vai\iF{g) for g £ £^ such that n[g] = n, (4-1) 
where F is defined by (12. Sp . is attained for a unique density operator g[n]. 
Proof. We denote 

A = {g £ such that n[g] = n} . 

Step 1: A is not empty. We start with a simple, but fundamental remark: thanks 
to our assumption on the density 77,(2;), the set A is not empty. Indeed, let (pi := 
||n||^y^-y/n and complete cj)i to an orthonormal basis (0i)i>i of 1/^(0, 1). The func- 
tion n belongs to H^g^. Hence, by Sobolev embedding in dimension one, n is con- 
tinuous and, from n > 0, we deduce that 

n(x) > minn > 
[0,1] 

and then y/n G Hpg^. For all ijj G L^(0, 1), consider the density operator v defined 

by 

i^i^:=V^{V^,i^), (4.2) 

we find 

i>l 

Jo 

Jo 

so, by the characterization (j2.2p . belongs to A. 



d 

dx 



dx < 00, 
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Step 2: F is bounded from below on A. From (|2.7p . we deduce that, for all g £ A, 
F{g) > n{x) log n{x)dx + ( ^^^^^'^^ ~ / "'(^)^^ > (4-3) 

since by Sobolev embedding n is bounded. Therefore, one can consider a minimizing 
sequence {Qk)k&i for (|4.ip . i.e. a sequence £ A such that 

lim F{Qk) = inf F(o-) > -oo. 

fc— ^+oo cr€^ 

5.' uniform bound in £. Let us prove that {okjk&i is a bounded sequence of 
Since E we already have 



\Qk\\ji 



Tr£>A:= / n{x)dx < +00. 
Jo 



Moreover, since the density operator defined by ()4.2p belongs to we have, for 
k large enough. 



Tr(a log - + TriVHQkVH) = F{gk) < F{v) + 1< +oo. 
Hence, using the inequality (j3.7p . we obtain 

-C7 {ti y/Hgky/H^ + Tr(^/;H'a/ff) < F{u) + 1 < +oo, 

thus 

supTr(\/^i?A;V^) < +00. 

fceN 

Step 4-' convergence to the minimizer. Since {gk)keN is a bounded sequence of 
one can apply Lemma [3. II to deduce that, after extraction of a subsequence, we have 

0k ^ Q hi i7i as /c — )■ +00 (4.4) 

and 

TT{^/Hg^/H) <limmfT^{^/HgkVH). (4.5) 
Next, by (j4.4p and Lemma [3.31 (ii), we get 

Tr{gk log - Qk) Tr(£)log £i - £») as +oo, 
which yields, with (j4.5p . 

F(£.) < liminf F(a) = inf F(cr). (4.6) 

Let $ E L°°(0, 1) and denote also by $ the bounded multiplication operator by <I>. 
Since gk converges to ^ in ^i, we have 

/ ^{x)n{x)dx = Tr{^gk) — )• Tr(<I>^) as k ^ +oo, 
Jo 

thus, from the characterization l\2.2\) . we deduce that n[g] = n, which means that 
g G A. This enables finally to conclude from (|4.6p that, in fact, we have the equality 

F{g) = inf F{a) = minF(cr). 

The uniqueness of the minimizer is a consequence of the strict convexity of F, see 
Item ( Hi ) of Lemma 13.31 □ 
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5. Characterization of the minimizer 

This section is devoted to the second part of our main Theorem I2.H the charac- 
terization of the minimizer. As we explained at the end of Section [21 we need to 
define a penahzed version of our minimization problem. 

5.1. A penalized minimization problem. Consider a density n{x) such that 
n > on [0, 1] and n E Hpg^. For all e G (0, 1] we define the penalized free energy 
functional, for all ^) S 



r7l|2 



Proposition 5.1. Let e G (0, 1) and let n G H^g^ such that n > on [0, 1]. The 

minimization problem without constraint 

min Fir (g) (5.1) 

where is defined above, is attained for a unique density operator [n] , which has 
the following characterization: we have 

g,[n]=ew{-{H + A,)). (5.2) 

where £ ^ler- 

Proof. Since the entropy functional Tr{glog g — g) is not differentiable on we 
regularize it. For all rj G [0, 1] and s G M^, we define the regularized entropy 

(3r,{s) = (s + ??)log(s + r/) - s-rjlogrj, 

and the associated free energy functional, for all g £ 

= Tr (/3^(£>)) + Tri^gVH) + ^\\n[g] - n\\l,. 

Notice that /3^(s) = log(s + rj), /3rf{0) = 0, and that is strictly convex on and 
liolomorphic on (— oo) x M for the convenient branch. 

Step 1: minimization of F^^^j- In this step, we prove that for all r] G [0,1], the 
problem 

min Fs^r^ig) (5.3) 
e&£+ 

admits a unique minimizer fj^^^. Notice that for r] = 0, this problem is nothing but 
()5.ip : in the statement of the Proposition, we have denoted shortly gs = gs,o- 
By (|2.6p and a Sobolev embedding in dimension one, we have 

||n[f?]||Loo <CTig + CTT{VHgVH), 

so the functional F^ is well-defined on for all r/ G [0, 1] and e G (0, 1]. We will 
need the following technical lemma on the function /3^. 

Lemma 5.2. The application g i— /3^(£») possesses the following properties. 

(i) There exists a constant C > such that, for all ^? G and for all r] G [0, 1], we 
have 

Tr/3^(^) > -C (TrVHgVHy^'^ . (5.4) 

(ii) Let g^ be a hounded sequence of £^ such that gj. converges to g in J\, then for 
all rj G [0, 1], P'q{gk) converges to I3^q{g) in J\. 

(Hi) For all t] G [0, 1], the application g i— )• Tr/3^(^j) is strictly convex on E^. 
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(iv) Consider a sequence q^i hounded in iS+ such that Qr^ ^ Q in Ji as rj ^ 0. Then 
Tr/3^(^^) converges to Tr/3o(i?) as — )• 0. 

Proof of the lemma. It is not difficult to adapt tfie proof of Lemma 13.31 in order to 
sliow Items (i), (ii), (Hi). We sliall only prove Item (iv), proceeding similarly to 
Step 2 of Lemma [3. 3 1 We first notice that the function /3.^ converges to Pq uniformly 
on all [0, M], M > 0, and that one has 

VsG[0,M], \P^{s)\<Cm^s, 

with Cm independent of rj. Let M = sup^ ||i?»)||£(L2) < +oo. For all N € N*, by 
using the inequality (j3.8p . we get 

\M\{q,])\ < Cm E y^S^ ^ {TT^e,^)'" < 

p>N p>N ^ ^ 

where we used the fact that {gk) is a bounded sequence of and where Ap[£»,^] 
denotes the p-th nonzero eigenvalue of g. Hence, decomposing 

p>N p>N 
one deduces from the uniform convergence of to /3o and from 

WIqt,] - >^p[q]\ < Wqti - e\\c{L'^) < Wqt, - qWji, Vp > i 

(see the proof of Lemma IA.2I in the Appendix) that 

\Tvpr,{Qv) -T^Me)] ^ as r/ ^ 0. 

The proof of Lemma 15.21 is complete. □ 
Let us now study the minimization problem ()5.3p . for fixed rj G [0,1]. For all 
g € we deduce from ()5.4p that 



Fe,^{g) > -C (TTVHgVHy^\TT{VHgVH) + U\n[g]\\l2 - ^ 



where we used (a — 6)^ > — 6^. From this inequality, we deduce two facts. First, 
that infgg£-_|_ Fg^^riio) > — oo. Second, that any minimizing sequence g^ is bounded 
in iS-L. Indeed, we have 



thus 



-C (riVHgkVHy^^ + TiiVHgkVH) + ^ sup ||n[a]|lii 

< supFe,^{gk) + Trll^llia ^ -|I'^IIl2- 

k Z£ £ 

This implies that, for all k, 

Ti gk + Tr ^/H gkVH < Ce (5.5) 
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where is a positive constant independent of k and rj. Hence, according to Lemma 
13.1^ one can extract a subsequence still denoted Qk such that the convergences 

Qk^ Q in Ji, Ty:{^^HQ^/H) < lim mi Tr {^/H QkVH) (5.6) 

fe— >+oo 

hold true as A; — )• +oo. By Lemma 15.21 (ii), we have 

Trf3^{Qk)^Tr/3^{Q). (5.7) 

Let us now prove that n[Qk] converges to n[g] in L°°(0, 1), which implies in par- 
ticular that 

W^iek] - n\\l2 \\n[Q\ - nllla. (5.8) 

We have 

2 

< TiQk + TrVHQkVH < +00, 



WVnigklWu^ = / n[gk]ix)dx + 



d 



L2 



where we used ()2.6p . Therefore, the sequence {^n[Qk\)k&i is bounded in H^(0, 1), 
and by Sobolev embedding one can extract a subsequence such that \/n[Qk\ con- 
verges to a function / E C^([0, 1]) in the L°°(0, 1) topology. This implies that n[Qk\ 
converges to in L°° . 

Moreover, for all <I> € L°°(0, 1), we know from (|5.6p that 



/ n[Qk\^dx = Tr q^^^ -^Tv = 
Jo Jo 



nlol^dx, 



which means that n[gk] converges weakly to n[g] in L^{0, 1). This enables to identify 
the limit: we have in fact /-^ = n[g]. 
Finally, 1^ and 1^ yield 

^e,»7(£') < liminf i^£,r?(^fc) = inf -Fe,r,(o"), 

so g G is a minimizer of ()5.3p . Furthermore, one remarks that the application 
^ I— )■ is linear, so the application 

Q i-> \\n[gk] - n\\l2 

is convex, and it can be deduced from Lemma 15.21 (in) that is strictly convex: 
the minimizer Qs,r][n] is unique. In the sequel of this proof, n being fixed, we 
denote shortly ^3^^^ instead of ^^^^[ri]. Notice that, from (|5.5p . one gets an estimate 
independent of the parameter 77 S [0, 1]: for all e £ (0, 1], one has 

sup Ti Qs^ri + sup Tr \/^i?£,»)V^ < +00. (5.9) 

r;G[0,ll ' velO,l] 

Step 2: differentiation of for r/ > 0. We will use the following lemma, whose 
proof is in the Appendix. 

Lemma 5.3. Let r] S (0,1]. Let g £ £^ and let cv be a trace-class self-adjoint 
operator. Then, the Gateaux derivative of the application 

g^F^{g) = TTf3r,{Q) 
at g in the direction uj is well-defined and we have 

DF^{g){u)=Tr{(3'^{g)uj). 
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Let h he a, bounded Hermitian operator. For g G consider the operator 
g + t^fgh^fg. Assume that t € [— toi^o]) with < toll^ll ^ 1- For such values of t 
and for all G L^(0, 1), we have 



> ll\/^'/'iii2(j^)( 

> 0. 



1 - Itl 



Therefore g + t^fgh^fg is nonnegative, self-adjoint and belongs to Ej^ since 
Tr{VH^h^VH) < = \\h\\TT{VHgVH) < 00. 



Moreover, we have the following estimates: 



\\V0h\ 



J2 



< 



\S\\ji, 



H^ghW 



J2 



< 



Ti{VHgVH), 



TTi^\h\^) < \\h\\\\g\\j„ TT{^^g\h\^g^) < \\h\\TT{^gVH). 



Therefore, by linearity, the following equality holds in Wp^{Q, 1) C L°°(0, 1) C 
L2(0,l): 



n[g + t^h^ 
which yields, for all t ^ 0; 
\\n[g + t^h^] - n 



n[g\ + tn[^h^, 



|2 

Il2 



n[^/gh^/g\{x) {n[g\ — n) {x)dx+0{t). 

(5.10) 

From Lemma 15.31 and from (IS.lOp . one deduces the following expression for the 
Gateaux derivative of F^^^ in the direction uj = ^/g^h^J'g^: 

Fe,r,{Qe,ri + t^/g^h^/g^^) - Fe^r,{Qe,ri) 



2t 



lim 



Tr {(3'^{ge,'n)^/g^h,/g^) + Ti {\fH ,/g^hy/gZ^\fH) + Tr(^e,^^/^/i^ 
Tr {^/g^ {/3'rj{Qs,n) + H + ^/g^h) , 



where we used the cyclicity of the trace and where we have denoted 

1 



(5.11) 



Ae^^{x) = -{n[ge,^]-n){x). (5.12) 

denotes here, with an abuse of notation, either the L°° function 
yle^^, or the operator of multiplication by Af,^^ which is a bounded operator. Indeed, 
A^^j^ belongs to Hpg^ C L°° since n G H^g^ according to the hypotheses and since 



Note that A^ 



£,r] 



Now, we have the tools to conclude: since g^^^^ is the minimizer of (|5.3p and since 
Q£,r] + ge^rjh^ ge,r] belongs to for t small enough, the Gateaux derivative (|5.1ip 
vanishes and for all h G C{L'^), self-adjoint, for all r] £ (0, 1], we have 

Tr (/3;(&,„) + H + A,,r,) Vs^h) = 0. (5.13) 



Step 3: convergence of g^^^i cls 0. From the estimate (|5.9p and from Lemma 
3.11 one deduces that there exists g G £^ (dependent of e) such that, as 77 — )■ 0+, 

ge,r, -^ginJi and Ty:{VHgVH) < liminf Tr(\/^e'e,r,\/^)- (5.14) 
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Then, from Lemma [5.21 (iv), one deduces that 

(&,,,) = l3o{g). (5.15) 

r/->-0 

Moreover, one can deduce from ()5.14p and from Sobolev embeddings in dimension 
one, exactly as to prove ()5.8p . that 



W^iQrj] - n\\l2 \\n[g\ - n||^2 
as ?7 — 7- 0. Together with (|5.14p and ()5.15p . this leads to 

F,,o{q) < lim Fe,r,{Qe,r,). (5.16) 
Moreover, by definition of q^^q and Qs^-q as minimizers of F(.^o and ^e,r], one has 

FefliQefl) < Feflio) and Fe,r,{Qe,v) < -^e,r;(&,o)- 
Applying Lemma 15.21 (iv), one gets 

lim Fe^niQe.o) = ^e,o(&.o), 

r)-)-0+ 

and finally all these limits are equal, since 

FefliQefl) < Fefiio) < lim Fe^r,{Qe,-n) < 1™ Fe^r,{Qefl) = -F'e,o(&,o)- 
r?-+0+ r)-i>0+ 

Hence, by uniqueness of the minimizer, we have q = g^^. Moreover, we deduce also 
from FsfiiQefl) = lim,,^0+-^£,7?(£'e,r;) that 

TiiVHg.oVW) = lim Tt{Vh nVH). 

Hence, by applying the second part of Lemma l3.1l we get finally 

Qe,rj — ^ ^,0 in Ji and VH^/g^ — > \fH^/g^ in J2 as — > 0. (5.17) 

Now we have the tools to pass to the limit in (|5.13p as — t- 0+. First, let us 
prove that, for all bounded operator /i, 

{(3'^{ge,ri)Qe,rih) = Tr (/3o(£'e,o)£'e,o ^) • (5.18) 
To this aim, we introduce a parameter k > and decompose 

Tr(/3^ (£»£,,,) &,r,/i) =Tr(lg^_ ,,>«;/3^(£>e,^)£>e,^/l)+Tr(lg^ ^<K/3^(£'e,^)&,r,^) . (5.19) 

Since /3^(s) = log(s + ??), on all interval [k, M] with < k < M , one has 

lim max |/3' (s) - /3^,(s)| = and \sl3'{s)-sl3'Js)\<Cs. 
v-^o se[K.,M] ' ' 

The first term in this decomposition (|5.19p can thus be uniformly approximated, 
for M large enough: 

|Tr {lg,^^>K /3'ri{Qe,v)Qe,v - Tr (lfe,o>K /?o(&,o)£'e,o h) I 

< ll&,»?ll^ill^ll max 1/3' (s) - /3o(s)| + C^Hfe,,, - eie.oll ll^ll, 

and converges to as — )• (for all fixed k > 0). Consider now the second term 
in the right-hand side of (|5.19p . We have a uniform bound s^/^/3^(s) < R for 



s £ (0,M], so 



^7/8 1 



Jl 
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1/4 



where we used again the bound ()3.8p for the eigenvalues Ap of ^j^^^, together with 
the estimate ()5.9p . Hence 

lini sup |Tr (lg^^<«,/3'(&,^)^e,^/i)| = 0. 

'^~*''Jr;G(0,l] 

This ends the proof of (jS.lSp . 

Second, by (|2.6p and by Sobolev embedding, we have 



1 



1 



{n[ge,o] - n) (x) 



in the L°°(0, 1) topology. Hence, the corresponding multiplication operators satisfy 



A, 



A,,o in£(L2(0,l)) 



and the convergence of jQe,r) in ^2 yields 



lim Tr (^/^A^,^ y/ol^h) = Tr [y/g^Aefl ^/g^h) . 



Third, the convergence of y H^fg^ in J2 yields 



lim Tr 



VH^^,h) = Tr (^VhVh 



(5.20) 



(5.21) 



Finally, one can pass to the limit in (|5.13p and (|5.18p . ()5.20p . (|5.2ip give, for all 
h S C{LP') self-adjoint, 

TV ( V& {\og{Qe) + H + A,)^eh) = 0. 

where we have denoted = g^^Q and A^ = Ai^fi. This means that 

^e{logige) + H + Ae)^e = 0. (5.22) 

Step 4: the kernel of g^ is {0}. In this step, we will prove that, for all e G (0, 1], the 
kernel of the minimizer of is {0}. 

Let us prove this result by contradiction. Assume that the kernel of g^^ is not 
{0} and pick a basis function (p G Ker^j^. We first complete <j) into an orthonormal 
basis {0, (^p)pg/} of Ker^£ (/ may be empty, finite or infinite). Then, we denote 
by (Ap)i<p<7v the nonincreasing sequence of nonzero eigenvalues of g^ (here N is 
finite or not), associated to the orthonormal family of eigenfunctions {(l)p)i<p<N- 
We thus obtain a Hilbert basis {(j), (ipp)p^f, {(j)p)i<p<N} of L^{0, 1). Since it is not 
clear whether <p belongs to H^g^, let us regularize it by setting 

= {l + aVH)-^(l), 

where a > is a small parameter. We have (p" G H^g^ and, as in the proof of 
Lemma 13.11 

lim 6°" ^ d> 
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in L2(0, 1). We simply fix a > such that |((/>°,0)| > 1/2. Denote by P° the 
orthogonal projection 

and consider the positive operator Q(t) = Qe + tP"' for t > 0. From (/>° G H^g^, we 
deduce that the operator belongs to We shall prove that there exists t > 
such that 

FMt)) < FeiQe), (5.23) 

which is a contradiction. 

Let ?7 > and denote as before /3(s) = s log s — s and ,S,,(s) = (s + r/) log(s + ??) — 
s — r/log ry. From the min-max principle and from the positivity of the operator P", 
one deduces that 

Vp G N*, \p{Q{i)) > Ap(&), 
where Ap(-) denotes the p-th eigenvalue of the operator. Hence, we have 

/3{XMt)) - = / log{s)ds 

< / log(s + 'r])ds 

which implies 

Ttimt)) - Tr(/3(&)) < Ti{(3r,{Q{t)) - Tr(/3^(^,,)) 

and then 

Therefore, to prove (|5.23p . it suffices to find rj > and t > such that 

Fe,Mt)) < FsAQs). (5.24) 
Since P" belongs to £^ and by Lemma |5.3| for all r/ > one can differentiate F^^^i 
at in the direction P" and one has 



t 



N 



One has 

Tr(log(£.,+7?)P") = W,ct>)\^\ogri + Y,WAv)?^og{Xp + ri) 

p=i 

+ ^|(<A°,Vp)l'log??, 
pel 

hence, by using |(i;^",i;^)| > 1/2, one obtains for < r] < 1/2 

1 ^° I 

TV (log(^, + r?)P") < - logr/ + 5^ , 0p)|2 log(Ap + -), 

p=i 

where po has been chosen such that Ap < 1/2 for p > pQ. Therefore, there exists a 
constant C^^a independent of rj (but depending on e and a) such that 

F.,,(g(t)) - F,,^(fe) ^ 1 , 

Im ^ '- < - log ?7 + Ce,a- 
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To conclude, it suffices to choose rj small enough such that j log rj + Cg^a < 0. Then 
one has 

t->0 t 

and for t small enough one has (I5.24p . which leads to a contradiction. This ends 
the proof of the claim. 

Step 5: identification of g^. Notice that, since A,. S H^g^ C L°°(0, 1), the operator 
Ha '■= H + Ag with domain D{H) is bounded from below and has a compact 
resolvent. Denote by (Ap, (/>p)pgN* the eigenvalues and eigenfunctions of g^. From 
the previous step, we know that, for all p, we have > 0. Moreover, {(pp)p^^* is 
a Hilbert basis of L^(0, 1). We will prove in this step that (c/ippgN* is a complete 
family of eigenfunctions of associated to the eigenvalues — logAp. 
Apply (|5.22p to (pp. Since from Step 4 we know that Ap > 0, we obtain 

^ (log(A^) + H + A,)cP'p = 0. 

Remark that, since \fH^fgl is bounded (with adjoint operator ^J^^/H), we know 
that (/)p belongs to Hpg,,. Taking the scalar product of the above equation with 



h^^ leads to 



(log(A;) + H + A,) 



= yA|log(A^)5p, + {yH(t>l,sfH^e4>i) + [As (p'p, VFefPl) 
= ^Al (log(A^)5p, + (/H^cA;, Vh^PI) + {Ae cp;, <Pl)) . 
Hence, for all p,q G'N*, 

^(Pl, Vh<PI) + {A, (PI, <Pl) = - \og{\l)5pq. (5.25) 



The family ((/ippg^* is thus an orthogonal family for the following sesquilinear form 
associated to Ha'- 

Qa{u, v) = (VHu, VHv) + {As u, v) . 



Note that, since A^ S L°°, there exists two constants m, M > such that 
Ii 1. 



VuGH^ , -^M^r <QA{u,u) + m\\u\\l2 <M\\u\\^^. (5.26) 



Let us now prove that this family {(pp)p,^fq* is dense in Hpg^. Let (p G Hpg,.. We 
already know that the following series: 



N 

P 

p=l 



(pN = ^{<t>l, 



converges in -L^(0, 1) io cp as N ^ +oo. We will prove that in fact this series 
converges in which, by (j5.26p . is equivalent to saying that 

Qa{(P,(P)= lim Qa{(Pn:(Pn)- (5.27) 
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Again, the key argument of the proof wih be the fact that Qe is the minimizer of 
F^: for all t > 0, we have 

Q<F,{Qe+tP)-F,{Qe) (5.28) 

where P denotes the orthogonal projection on (j): 

Pu:= (l){(l),u), VnGL2(0,l). 

Indeed, (j) G H^g^ implies that P S thus Q{t) := + tP belongs to £^ for all 
t > 0. Now, as in the previous Step 4, one can prove that, for all r/ > 0, we have 

< F,{0, + tP) - F,{q,) < Fe,^{Qe + tP) - F,^r,{Qe) 

and 

^:^^ Fe,n{Qe + tP)-Fe,r,{Qe) ^ Tl {\og{Qe + ll)P) + Tv{VH P^) + Tt{A,P) 
t^O t 

= |(<A,'APl'log(A; + r/) + QA(</','A)- 

pen* 

Therefore, for all r/ > 0, one has 

- E l('A,0Pl'log(A; + r/)<QA(</',<A). (5.29) 

pgN* 

Let N &N* large enough, such that log(A^) < 0. For rj > small enough, one has 

Vp > N, log(A^ + ??) < 0, 

thus (fOOjl yields 

N 

-^\{4>,<^l)\^log{Xl + v) <Qa{^A)- 

p=l 

Since we know that A^ > for all p, one can pass to the limit in this inequality as 
77 ^ 0: 

N 

-Y,\{(t>,^l)\^logXl < Qa{^A)- 
p=i 

Remarking that, by the orthogonality property (j5.25p . one has 

N 

Qa{(^n,M = -^\{(t>,4>l)\^^ogXl, 
p=l 

this inequality reads 

Qa{^n,M <Qa{^,^). (5.30) 
In particular, this means that ((/>Ar) is a bounded sequence of H^g^, thus converges 
weakly to </) in as — )• +oo. From the equivalence of norms (j5.26p . we then 
deduce that 

Qa{4>,(I)) < liminf Q^((/)Ar,0Ar). 
A'^— >-+oo 

Together with (|5.30p . we get (j5.27p and our claim is proved: c/^at converges to (p in 
the H"*^ strong topology and the family {(pp)p,^^* is dense in H^g^. 

This enables to conclude the proof. Indeed, this density property implies that 
()5.25p is equivalent to 

V0 G H^g„ (Vh^;, yUA + {A, (p) = -{4>l, 0) log a;. (5.31) 
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This means that {(pp)p£f^* is a complete family of eigenfunctions of Ha (still identi- 
fied with the associated quadratic form) and that the associated eigenvalues of Ha 
are — log A^. In other words, we have 

Qe = exp {-Ha) 

in the sense of functional calculus. The proof of Proposition 15. II is complete. □ 

5.2. Passing to the limit. In this subsection, we terminate the proof of our main 
Theorem 12.11 Let such that n > on [0,1] and ^/n £ H^g^. By a Sobolev 

embedding in dimension one, -^/n is continuous on [0, 1], so that we have 

< m := min n(x). (5.32) 

xe[o,i] 

Propositions 14. l l and 15.1] define the unique minimizers Q[n] and Qeln] - shortly de- 
noted Q and here - of the minimization problems (j4.1l) and (|5.ip . Moreover, 
takes the form (|5.2p . with 

A, = ^{n'-n)eill^r, 

where we have denoted := Let us study successively the limits of g^, 

and as e — )• 0. 

Step 1: convergence of q^. Recall that, for all o" E we have 

F,(a)=F(a) + ^||nH-n||22 
and that n[^>] = n. Hence, by definition of q^^ we have 

< F,{q) = F{q). (5.33) 



Therefore, one deduces from the estimate (j3.7p that y H Qi^y H is bounded inde- 
pendently of e and that n[f)e] converges to n in -L^(0, 1). In particular, by Cauchy- 
Schwarz, we obtain 

-1 



|Tr — Tr q\ 



{n[Qi,] — n){x)dx 



< \\n[gs] — n\\i2 — )• as e — )• 0. 



The family is thus bounded in £^ independently of e and then, by Lemma 13. 1| 
there exists g € such that 

ge-^gin Ji and Tr{VHgVH) < limmf Tt{Vh g^VH). (5.34) 

Therefore, by Lemma 13.31 (ii), by the expression (|2.5p of F and by (|5.33p . one gets 

F(g) < liminf F(^s) < F(g). (5.35) 

Furthermore, we have n[g\ = n. Indeed, the strong J'l convergence of g^ implies the 
weak J'l convergence, thus 



V99GL~(0,1), [\[ge]{xMx)dx = Trig, Trig if) = C 

Jo Jo 

and we already know that n[gg] converges to n in L^(0, 1), so 



n['g\{x)(p{x)dx 



Vv? G L°°(0,1), / n[g\{x)^{x)dx= [ 
Jo Jo 



n{x)ip{x)dx 
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and then = n. 

Finally, ^ is a minimizer of (j4.ip and the uniqueness of this minimizer proved in 
Proposition 14.11 vields Q = Q. Moreover, one has F{Qf,) — )• F(q), so 



lim F{Qe) - limTr(£»elog & - Qe) 

e-5>0 e-s>0 



= F{0) -Tr{glogQ- g) =Tr{VHgVH). 
The second part of Lemma 13.11 can thus be applied and one has finally 
gs ^ g in J\ and ^fH^/gl — )• ^/H^ in Ji as e — >■ 0. 
Notice that, from Lemma [3.11 we also have 

\/& ~^ ^Tq ill Ji as e — )• 0. 



(5.36) 
(5.37) 



Ste'p 2: convergence of := n[iJ£]. In this step, we will prove that 

lim ling — nlljji = 0. (5.38) 

£-S>0 

From the previous step, we know that ge ^ Q in J\. According to Lemma |A.2| 
for the sequel of the proof, let us choose some eigenbasis (^ppeN* and ('/'p)pGN*j 
respectively of g^ and g and such that 



VpGN*, lim 



0. 



We claim that 



lim 



Qs(p% - \fQ(\>v\\ .o = 0, 



and lim ^fHJ^cff - 
e-s>0 ^ 

Indeed, to prove (|5.40p . it suffices to write 



H^/g4>p 



L2 



(5.39) 

(5.40) 
(5.41) 



L2 



+ 



< 
< 



'Qe\\c{L^) \\9p - 9p\\i^2 



^ - Ve)^p\\L^ 



'Qe\\j2 



+ 



J2 



h - ^/^|| 

then to use (|5.37p and (j5.39p . To prove (|5.4ip . one writes similarly 

^yfge¥p-^^fQ4>p „ < II^V^I|j2 ll'/'p -</'p||i2 + "/H y/gl - ^ y/g 



J2 



then use ([QGll and (f09]l . 

Let us prove (j5.38p . We already have — )• n in -L^(0, 1), so it remains to prove 
the convergence of in L^. One has 



n. 



pGN* 



thus 



driir dn 
dx dx 



< 



pgN* 

peN* 



^ dx 



^'p4>l 



pYp 



Therefore, from (j2.ip and Cauchy-Schwarz, one gets 
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with 





dn 


dx 


dx 







L2 



< sf + sl 



Hy/Q)(t)p 



L2 



peN* 



(5.42) 
(5.43) 



Let us use again ()3.8p . From the Gaghardo-Nirenberg inequahty, one deduces that, 
for all iV G N*, 



p>N p>N 

\p>N J 

1/4 / \ 1/4 



and similarly 



\P>N J \P>N^ J 



1/2 



< ^(Tr,,)V4(TrVF,,VF)'^'< ^ 



iVl/4' 



Ari/4^ 



E ^p\\^p\ 

p>N 



< 



c 



Thus, for ah iV G N* 



Ari/4- 



^pGN* 

c 



\\^^e<P'p-VH^4>p\\l2 

\p<N 



J2 (\\VHVp'e^;\\i2 + wVh^ci^pWI, 

p>N 

— ^ 1,2 C 

H^/p(t>p\\L2 + 



p<N 



iVl/4 



and from (|5.4ip one deduces that 



lim Isf 1^ = 0. 

£-S>0 



(5.44) 



26 



F. MEHATS AND O. PINAUD 



Furthermore, 

i^pGN* ) \p<N 

+ [ Y W^^QUI^ I ( E {\\VP'e(t>l\\l^ + \\V~P<t>v\\l^) 

\pm* ) \p>N 

ypm* / \p<N 

\ 1/2 

E II^^^/p7<^; - Vh^<Pp\\1. + E Il^^v^'/'plli2 



yp<N J pGN* 



X 



where we used again the Gaghardo-Nirenberg inequahty. Hence, from ()5.4ip . one 
deduces that 

hm = 0. (5.45) 
The convergence (j5.38p of is proved. 

Step 3: convergence of A^. By a Sobolev embedding in dimension one, Hpg^ is a 
Banach algebra: for all u,v £ Uper the product uv also belongs to H^g^. Hence, 
from (|5.3ip . one gets 

vp G r , V,/. G Hi,„ (VF<^^, VHi<p;cj,)) + = -{^;, 4>'p^) log a;. 

Multiply this identity by and sum up on p. Since = X^p-^pl'^pP' obtain 
that, for all (p £ H^^^, 



JO 



A,n,(j)dx = - E {^p,<PiQe^OgQe)(l)p) - E (^V&<^p) V^(</'\/&</'p)) 

pSN* peN* 

= - Tr(0(& log &)) - E [VH^e'e^'p, VH{ct)^,ct>l)) . (5.46) 

pGN* 

Let us examinate separately the convergence of the two terms in the right hand-side. 
From (|5.36p . Lemma [3.31 (ii), and from 

\TT{(t){Qelogge)) -Tr{(t){Qlogg))\ < \\ct>\\L°-\\gelogge - gloggWj^ 

one has 

lim sup \Tr{(l){gelogge)) -Ti{4>{g\ogg))\ =0. (5.47) 

^^°II<^IIhi<i 

Let us now prove that 

^(PVOI VH(P^ in J2 as e ^ 0, (5.48) 
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where (j) denotes the operator of muhipUcation by </). Using the identification (|2.1 
one gets 



E 



2 

L2 



< 2 



< C 



dx 



L2 



El 

p 

E 



L2 



L2 



H{^/g^ - ^/g)(t>p 



L2 



2 

L2 



where we used a Gaghardo-Nirenberg inequahty. Hence, from (|5.36p and (|5.37p . one 
deduces ()5.48p . Finahy, from ()5.36p and (|5.48p . one gets the fohowing convergence, 
as e — )• 0: 



J2 



J2 



(5.49) 



pGN* 

Let us now define a linear form on Hpg^. For if) G H^g^, we set 

(^'V')h-,hi.. = Tr (^(^'log^.)) + E (^^/^'Ap)) . (5.50) 



From the above estimates, one deduces that 

n jji 



< C 



Tr log q\ + Tr + Tr \/^ qVh) ■ 



(5.51) 



Since n(x) > m > on [0, 1], the application ip ^ ^ is continuous on H^g^,, so the 
above defined linear form A is continuous on H3,g^ and belongs to its dual space 



Hpg|,. Moreover, we have proved by ([06]) . (fOTP and (f5^ that, for all G 



pen 



lim sup 



e-5>0 



0. 



To conclude, it remains to use the convergence ()5.38p of to n in H^, which implies 
that — converges to - and that, in fact. 



lim sup 



Aetpdx - {A,^)^~i 
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In other words, one has 

^ A in the H~g^ strong topology as e — t- 0. 



(5.52) 



Step 4-' identification of g and conclusion. Let us define the following forms on H^g^, 
QAA^,ij) = {VHip,VHi;) + {A,ip,il;) 



The form [A, (pip) symmetric, form-bounded perturbation of ( v H(p, v Hip) 

with relative bound < 1. Indeed, by (|5.5ip and by ^ G H,^ 



per' 



(AI^/^Hh-ihi ^ C|||<^|2||hi <C||v.||i4 + C||(/.||i. 



dip 



dip 



dx 



L2 



+ C\\ip 



1 1/2 



dip 



dx 



dx 

3/2 
L2 



L2 



< ^(^/F(^,^/:H^(/,) + C||(^||i2, 



where used the Young inequality and a standard Gagliardo-Nirenberg inequality. 
Let Ha be the unique self-adjoint operator associated to Qa- Then, according to 
Theorem XIII. 68 of |31) Ha has a compact resolvent and we denote by (/ip)pgN* its 
eigenvalues. In addition, we have 

\QAA^,V^)-QAi^,V^)\ < 



lHi> 



< C\\A, - A\\^-i ({VHip, y/Hip) + M\l2) . 



Moreover, Theorem 3.6 of [25j, chapter VI, section 3.2 yields the convergence of 
operators in the generalized sense, which implies in particular the convergence of 
the eigenvalues: 

- log a; = ^ ^ip yp G N* 
as e — )• 0. Hence, by continuity of the exponential function, 

Xp = exp{-fip) — > exp(-^p) Vp G N*. 

Besides, according to Lemma IA.2I in the Appendix, the jTi convergence of to g 
implies the convergence of the eigenvalues: 



Xp 



Xp yp e N*. 



This enables to completely identify the eigenvalues of g: we have 

Xp = exp(-/ip). 

Furthermore, from (I5.4ip and from Xp > 0, one deduces that (j)p — )• (j)p 
can thus pass to the limit in (|5.3ip : for all p G N* and for all (p G H^g^,, 



in H^. One 



which yields 

Qa{4>p, (p) = fj-pifpp, fp), ^4> e Hpg^, 

so, finally, (</>p)pgN* is the complete basis of eigenvalues of Qa- We have completely 
identified g: 

g = eicp{-{H + A)) . 



The proof of our main Theorem 12.11 is complete. 
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Appendix A. 

Lemma A.l. Let q £ and denote by {pp)p>i the nonincreasing sequence of 
nonzero eigenvalues of q, associated to the orthonormal family of eigenf unctions 
{(j)p)p>i. Denote by {Xp[H])p>i the nondecreasing sequence of eigenvalues of the 
Hamiltonian H . Then we have 

Tt{^pVh) = Y,Pp {Vh^p, VHcPp) >Y,Pp \p[H]. 

p>i p>i 

Proof. Notice first that 

p>i 

= '^Pp{VH(t)p,\^(t>p)- 

p>i 

Then, 

N N N-l 

Pp {VH(pp, Vll(j)p) = PN ^(V^0p, VHcpp) + (pN-i - Pn) ^ {\^4>p, "/Hfpp) 

p=l p=l i=l 

2 

+ • • • (P2 - P3) VH(t>p) + PliVW^l, VH^i). 

p=l 

Using [28) . Theorem 12.1 page 300, it comes 

N N 

p=i p=i 
so that, since pp < Pp^i, Vp > 1, 

N N N-l 

Pp {y/H(j)p, \fH(t)p) > PN ^p[H] + {pN-1 - pn) Ap[jj"] 

p=l p=l p=l 

2 

+ ---{p2-pz)Y.^p[H]+piXi[Hl 

p=i 

N 

p=l 

We conclude by passing to the limit as — t- +oo. Notice that the theorem of |28) 
is written for Hamiltonians defined on M.'^, but it can be easily extended to bounded 
domains. This ends the proof of the lemma. □ 



Lemma A. 2. Let a sequence converging to q in Ji as k ^ +oo. Then the 
corresponding nonincreasing sequence of eigenvalues (Ap)pgN*; {\)p&N* converge as 
follows: 

yp G N*, lim ^ Xp. 

fc— >-+oo 
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Moreover, there exist a sequence of orthonormal eigenbasis ((/>p)peN* of and an 
orthonormal eigenbasis {(j)p)p^^* of g such that 

VpGN*, lim U'; -<Pp\\ L2=0. 

Proof. Let us first prove the convergence of the eigenvalues. According to |33| . 
Theorem 1.20, we have the fohowing relation between the eigenvalues of gf^ and g: 



oo 

Ap - Ap = ^ apgXg[gk - g], P>1, 

9=1 



where {Xq[gk — Q])q>i denote the eigenvalues of gk — Q and a is a doubly stochastic 
matrix, that is a matrix with positive entries such that YlpLi ^pq — Tl'qLi ^pq — 1- 
The minmax principle [3T] implies |Ag[£)fc — ^]| < \\gk — q\\c{l'^)i W ^ 1; so that 

|Ap - Apl < \\gk - q\\c{l^) < \\gk - qWji, p>'^, 

which gives the desired convergence property. 

Let us now prove the convergence of eigenfunctions by following |25) and |2]. Let 
cr(g) be the spectrum of g and consider an eigenvalue Ap of p with multiplicity rup. 
Let dp be the distance between Ap and the closest different eigenvalue, 

dp = min |Ap — iJ,\, 

and denote by T the circle of radius ^ centered at Ap. Assume k is large enough so 
that 

\\Qk-Q\\c{L'^)<Y- 

Then according to [25], theorem IV. 3. 18, (see also example 3.20), there are exactly 
TTip (repeated) eigenvalues of gk included in F and denote by 0^ / = 1, • • • , mp the 
associated eigenfunctions. If {(t>p,i)i=i,--- ,mp denote the eigenfunctions associated to 
Ap, we construct an operator Up such that 

4>p,l = Up4>p,i, 1 < / < rrip, and Up —?■ I in C{L^) as A; oo. 

To do so, let Pp[g\ be the projection operator onto the spectral components of g 
inside F, 

1 



According to [25], II. 4. 2, Remark 4.4, if 

\\Pp[Q\-Pp[Qk]\\c{L^)<l, (A.2) 
then an expression of Up can be given by 

U^ = (/ - (Pp[a] - Pp[Q]?y^' {Pp[Qk]Pp[Q] + Pp[Qk]){I - Pp[Q])) . (A.3) 
Let us verify first that ()A.2p holds for k large enough. We have 

Pp[Q]-Pp[gk] = ^ j^iizl - g)-' - {zl - gkr')dz, 

= / {zl - gy^{g - gk){zl - gky^dz. 
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From the definition of dp, we have 



sup \\{zl - g)-'\\^L^) = —. (A.4) 
zer dp 

Moreover, owing (lA.ip and noticing that 

{Zl - Qkr^ = {I + {Zl - 0)-\q - Qk))-\zl - q)-^ 

and 

\\{I+{zI-e)-\Q-Qk))-^\\^^^^.^ < {i-\\{zi-e)-\Qk-Q)\\aL^)y\ 

< (^^- ^\\Qk- eWciL^)^ 

we conclude that 



-1 



SUp\\{zI-Qk) '\\ciL^)<^ II ^ II ■ 

^er -f -\\8k- gWciL^) 

Therefore, (|A.4p and (|A.5p imply the inequality 



(A.5) 



\Pp[g] - Pp[Qk]\\c{L^) < 



f - WQk - q\\c(l2) 



Assuming k is large enough so that \\Qk — q\\c{l'^) < ^) the above inequality yields 
the desired result since 



\Pp[Q] - Pp[Qk]\\c{L^) < -r\\Qk - q\\c(l2) < 1- 



(A.6) 



Let us prove now that Up — t- / in C{L'^) as A; — t- cxd. First, remarking that Pplg] = 
Pplol'^ and -Pp[gfc] = Pplgk]"^ since both are projections, ()A.3P can be recast as 

= {l- (PpiQk] - Pp[g]fy~' {I + Pp[gk]iPp[Q] - Pp[gk]) + {Pp[gk] - Pp[g])Pp[g]) • 

Let 5 := \\Pp[Q] - Pp[gk]\\c{L2) < 1- Then 
and, together with ()A.6p . 



< (1 - 5' 



< (1 - (\\Pp[0k]{Pp[g] - Pp[gk])\\c(L^) 



+ WiPpiQk] - Pp[g])Pp[Q]\\ciL^) + WiPpiQk] - Pp[g]y\\ciL^)^ 
< Cp\\gk - g\\c{L2), 

where the constant Cp does not depend on k for k large enough. It thus follows that 
\\(f>p,i - 4>p,i\\l^ ^ CpWQk - g\\c(L2) ^ as A; ^ oo. 



This ends the proof of the lemma. 



n 



^.riig + tOj) = ^ I Py{z){zl - Q-tu) ^dz. 
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Proof of Lemma \5.3[ Let T be an oriented curve in (— ??, +00) x R that contains the 
interval (— ^,2||^||). Let t G [— toj^o]) with 2to||w|| < min(r/, ||f)||). For such values 
of t, the spectrum oi g + tu is included in the interval (— ^,2||^||). 

Since /3,j is holomorphic in (— ry, 00) x R, one can define /3^(£) + toj) and f3rf{g) by 

/3r,(£') = ^ I'^f3r,iz){zl - g)~'^dz, 

1 
2in 

Let |t| E [0, min(to, ii)], where 

ii dist(r,a(£>))"^ ll^ll < 1, 
(T{g) denoting the spectrum of g. We have 

{zl - g - tuj)-'^ - {zl - g)-^ = {zl - g)'^ {I - tuj{zl - g)-^y^ - {zl - g)-^ 

km* 

where the latter serie is normally converging in Ji. Indeed, first, 
\\{tio{zI - gr'fWj, < \tnu:{zl - g)-%^, 
so that we only need to estimate ui{zl — g)^^ in Ji. Now, since g is self-adjoint, 
\\{zl - g)-'^\\ = d\si{z,a{g))-^ < dist{T, a{g))-^ , 

and it comes 

\t\MzI - g)-^\\j, < \t\\\{zl - g)-^\\Mj„ 

< \t\dist{T,a{g))-^\\uj\\j^, 

< 1, y\t\<ti. 

We thus can write 

t'^ {{zl - g- tcj)"^ - {zl - g)~^) = {zl - g)~^oj{zI - g)~^ + tA{t, z), 

where the operator A{t, z) is uniformly bounded in Ji with respect to t and z for 
\t\ G [0, min(to, ii)] and z G F. Hence, 

t-i[Tr/3^(^ + tu;)-Tr/3^(f?)] = j j3^{z){zl - g)-^uj{zl - g)-^dz 

+^Tr^/3^(z)^(t,z)dz. 
The two expressions of the right-hand side are well-defined and we have 
Tr / ^^{z){zl - gY^u){zl - gY^dz 



< \Pr,{z)\\\{zI - g)-^u{zl - g)-^\\j,dz, 

< distiT, a {g)y^\\uj\\ [ \P^iz)\dz, 



Tr / /3^{z)A{t,z)dz 



< sup \\A{t,z)\\ J, / \f3^iz)\dz < C\, 
zer 
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where Ci is independent of t E [0, min(to, ^i)]- Hence 

lim rMTr + - Tr/3„(£i)l = —Ti (3Jz)izI - gV^ujizI - g^^dz, 
t-)-o+ 2z7r Jy 

iGN* ^ 

where (pi, denote the spectral elements of g. Standard complex analysis then 

implies that 

^ f3vi^)iz - PiY'^dz = (3'^{pi). 
We therefore get the following expression of the Gateaux derivative: 

DF^{g){u) = J2 /5;(Pi)('Ai,^0i)- (A.7) 

iGN* 

The serie is absolutely converging since /3^(s) = log{s + rj) is locally bounded on ]R-|_ 
(recall that r/ > 0): 

\f3'^ipi)i<Pi,u;cl)i)\<C\icl)i,u;<Pi)\, 

and since we have assumed that G jTi. Finally, to identify the derivative, it suffices 
to notice that 

P'^{pi){(Pui^(t>i) = Y.{(Pi,u:fi'^{gm = Tv{u:fi'^{e)) = TT{(3'^{e)u). 

iGN* iGN* 

The proof of Lemma 15.31 is complete. □ 
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